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Quantum interfaces between photons and ensembles of atoms have emerged as powerful tools for 
quantum technologies. A major objective for such interfaces is high fidelity storage and retrieval of a 
photon in a collective quantum state of many atoms. This requires long-lived collective superposition 
states, which is typically achieved with immobilized atoms. Thermal atomic vapors, which present 
a simple and scalable resource, have, so far, only been used for continuous variable processing 
or for discrete variable processing on short time scales where atomic motion is negligible. We 
develop a theory based on the concept of motional averaging to enable room temperature discrete 
variable quantum memories and coherent single photon sources. We show that by choosing the 
interaction time so that atoms kept under spin protecting conditions can cross the light beam 
several times during the interaction combined with suitable spectral filtering, we erase the “which 
atom” information and obtain an efficient and homogenous coupling between all atoms and the 
light. Heralded single excitations can thus be created and stored as collective spinwaves, which 
can later be read out to produce coherent single photons in a scalable fashion. We demonstrate 
the feasibility of this approach to scalable quantum memories with a proof-of-principle experiment 
with room temperature atoms contained in microcells with spin protecting coating, placed inside 
an optical cavity. The experiment is performed at conditions corresponding to a few photons per 
pulse and clearly demonstrates a long coherence time of the forward scattered photons, which is the 
essential feature of the motional averaging. 

PACS numbers: 32.80.Qk, 42.50.Ex, 42.50.Pq, 03.67.Bg, 42.50.Pq, 03.67.Hk 


Quantum systems can potentially enable powerful 
quantum computation m and highly secure quantum 
networks [IHZl- Especially for the latter, it is essential 
that the information can be stored in quantum memories 
for processing [ill]. To this end ensembles of cold atoms 
have previously been considered for quantum memories 
since their large number of atoms enables a strong light- 
atom interaction IMH]. Cold atoms, however, require 
extended cooling apparatus, which makes the scalability 
of such systems challenging. In contrast, room tempera¬ 
ture atoms are much simpler to work with and have been 
used for a range of operations with continuous variables 
[S]. Several experiments with room temperature atomic 
cells have exploited a form of motional averaging where 
atoms move in and out of the beam many times during 
the interaction. By using a spin-preserving coating of the 
walls of the cell [5], the atoms can return back into the 
beam after colliding with a cell wall without losing the 
phase information. The averaging thereby removes the 
detrimental effect of atomic motion for continuous vari¬ 
able quantum information processing, where a specific 
optical mode is measured by homodyne detection m- 
m- On the other hand, for discrete variables protocols 
based on the detection of photon clicks, the situation is 
different since photon counters select an almost instanta¬ 
neous temporal mode. As a consequence, the efficiency 
of such systems for discrete variables where e.g. a sin¬ 
gle spin excitation is stored collectively in the ensemble 


is still limited by the incoherent atomic motion, which 
leaks “which atom” information and collapses the collec¬ 
tive state m- 

Here, we introduce room temperature microcells as a 
new system for discrete variable, ensemble-based quan¬ 
tum information processing. We show theoretically how 
the detrimental effect of atomic motion can be circum¬ 
vented in order to have an efficient and coherent inter¬ 
action between an atomic ensemble and light at the sin¬ 
gle photon level. The proposed technique can be used 
to make efficient single photon sources with memories 
and offers a solution to scalable photonic networks based 
on room temperature atoms. The spin-protecting micro¬ 
cells investigated here can also be used for long distance 
quantum communication in DLCZ-like repeater protocols 
[SIEIIIH] or quantum simulations [IMI]. The scalability 
of the system compared to cryogenic or cold atoms sys¬ 
tems opens up the possibility of employing large arrays 
of such systems combined with spatial multiplexing to 
enhance the communication rate [iiin]. As opposed to 
previous ensemble-based experiments with discrete vari¬ 
ables encoded in moving atoms |2J) , which typically rely 
on performing operations sufficiently fast that the atoms 
remain inside the laser beams, we show how a form of 
motional averaging similar to the one used for continous 
variable processing can be used to erase the “which atom” 
information. This novel approach to photon counting ex¬ 
periments alleviates the effect of atomic motion and can 
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FIG. 1. Sketch of (a) the atomic level structure and (b) the 
experimental setup for creating a symmetric Dicke state, (a) 
All atoms are initially pumped to state |0). The transition 
|0) —>■ |e) is driven by a weak laser held (fl) while the cavity 
mode (gr) couples |e) and |1). (b) The atomic ensemble is kept 
in a small cell inside a single-sided cavity with a low hnesse 
(cell-cavity). The quantum photons (thin arrows) are coupled 
from the cell-cavity into a high hnesse cavity (hlter-cavity), 
which separates them from the classical held (thick arrows) 
and averages over the atomic motion. Finally, the quantum 
photons are measured with a single photon detector (SPD). 


be seen as “trapping” coherent spins in a box-like po¬ 
tential made up of the coated walls. We also present a 
proof-of-principle experiment demonstrating the effect. 
Besides the specific experimental realization described 
here, the ideas behind this form of motional averaging are 
generally applicable and may be used in other systems 
where fluctuations of the coupling strength is an issue, 
such as ion crystals [24] . A related motional averaging of 
frequency fluctuations has previously been considered in 
super-conducting qubits |5S|. 

We consider a setup where an ensemble of atoms with a 
A-scheme level structure is kept in a small alkene coated 
cell Uni HZ] (see Fig. 10). A cell with quadratic cross 
section with side length 2L = 150 /im containing Cesium 
atoms was used in Ref. [5S|, for which the average time 
between atom-wall collisions was ^1.4 gtsec and the co¬ 
herence time was 10 ms. The atoms can thus endure 
several collision with the wall before loosing coherence 
making the cells suitable as quantum memories. The 
ensemble is kept at room temperature and, to enhance 
the interaction with the light, the cell is placed inside a 
single-sided optical cavity (ceZZ-cavity). In the proof-of- 
principle experiment (see below) a finesse oi T = YJ has 
been set by the output mirror transmission of 20% and 
the reflection losses on the cell windows but a cavity with 
a higher finesse F = 100 can easily be envisioned. The 
light leaving the cell-cavity is coupled into another high 
finesse cavity (/iZter-cavity), whose purpose is described 
below. 

Initially, all atoms are pumped to a stable ground 
state |0) (see Fig. 10). In the write process, the ob¬ 
jective is to create a single, collective excitation in the 
ensemble^ thereby creating the symmetric Dicke state 
|'0d) = S'd|00...0), with Sd = |l)i(0| where j 

is the atom number, N is the total number of atoms and 
|1) is another stable ground state in the atoms. The 
|0) —>■ |e) transition is driven with a laser pulse, which 


is far-detuned from the atomic transition to suppress the 
effect of the Doppler broadening of the atomic levels and 
absorption. In addition, the pulse should be sufficiently 
weak such that multiple excitations in the ensemble can 
be neglected. The write process is conditioned on de¬ 
tecting a single photon (quantum photon) emitted in a 
Raman transition |0) —)■ |e) —)■ |1). Upon detection, the 
atomic state is projected into the symmetric Dicke state if 
the light experienced a homogeneous interaction with all 
atoms, i.e. if the probability for different atoms to have 
emitted the photon is equal. In a realistic setup, the laser 
beam does not fill the entire cell and only atoms that are 
in the beam contribute to the cavity field, resulting in an 
asymmetric spin wave being created. Atoms leaving the 
beam will, however, return to the beam due to the fre¬ 
quent collisions with the cell walls. During the collisions, 
the atomic state is preserved due to the alkene coating of 
the cells and we exploit this to make a motional averaging 
of the light-atom interaction. If the interaction time is 
long enough to allow the atoms to move in and out of the 
beam several times, they will on average have experienced 
the same interaction with the light. Consequently, the 
detection of a cavity photon will, to a good approxima¬ 
tion, project the atomic state to a Dicke state. Since the 
cell-cavity has a limited finesse, it may, in practice, not 
have a sufficiently narrow linewidth to allow this averag¬ 
ing. We therefore introduce an external filter-cavity. As 
we show below, the output from the cell-cavity consists 
of a spectrally narrow coherent component and a broad 
incoherent component (see Fig. [^. By selecting out the 
coherent part, the filter-cavity effectively increases the in¬ 
teraction time and allows for motional averaging. At the 
same time, the filter-cavity can also separate the quan¬ 
tum photon from the classical drive if there is a small 
frequency difference between the two such that only one 
frequency is resonant in the filter-cavity while both are 
sustained in the cell-cavity. Furthermore, choosing the 
difference to be an even number of free spectral ranges of 
the cell-cavity ensures an overlap of the field modes at the 
center of the cavity such that atoms can interact simulta¬ 
neously with both modes if the length of the microcell is 
small compared to the wavelength corresponding to the 
frequency difference m- 


After a successful creation of an excitation in the en¬ 
semble, the state can be kept until it is read out. In 
the readout process a long classical pulse addresses the 
11) —>■ |e) transition such that the single excitation is con¬ 
verted into a photon on the |e) —>■ |0) transition (Fig. 10 r 
with g and D interchanged). This pulse should be long 
enough to allow for motional averaging as in the write 
process. The filter-cavity can once again be used to fil¬ 
ter the quantum photon from the classical drive photons. 
Furthermore, it can also be used to filter away incoherent 
photons as described below. 
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FIG. 2. Experimental data and simulation of the power spec¬ 
tral density (PSD) of the emitted light. The broad feature 
originates from short-time, incoherent light-atom interaction 
while the sharp peak is the long-time, coherent light-atom in¬ 
teraction. The optical scattering was obtained through the 
Faraday effect (see Sec|^for details) and is centered around 
the Larmor frequency at 823.8 kHz (the single high point in 
the hgure). The simulations have been rescaled to coincide 
with the data at this point. In the simulations we include the 
possibility of atoms being trapped in the coating of the cell 
walls. From the figure, we estimate that such trapping time 
is below 0.1 fj,s and can thus be ignored. The cell used in the 
experiment had dimensions 2L x 2L x 2Lz with L = 150 fim 
and Lz = 0.5 cm and the light beam had a Gaussian profile 
with a waist of 55 fim. The statistical uncertainty of each 
experimental point is very small and is therefore not shown. 


I. WRITE PROCESS 


To characterize the quality of our system, when con¬ 
sidered as a single photon source with memory, we first 
derive the efficiency of the write process and later discuss 
the readout efficiency and quality of the single photons 
being retrieved. The Hamiltonian, describing the write 
process, is {h= 1) 




i=i 


where A = wiaser — Weo with wiaser being the frequency of 
the driving laser and coeo being the transition frequency 
between the levels |e) and |0). We have assumed that 
the cavity is on resonance with photons emitted on the 
|e) —)■ |1) transition (see Fig. 10 1 ). flj (gj) character¬ 
izes the coupling between the laser (cavity) field and the 
j’th atom. The field in the cell-cavity is described by 
the annihilation operator Sceii and we have defined the 


atomic operators amh = \m)j{n\ for the j’th atom, where 
{m, n} S { 0 , 1 , e}. 

From the Hamiltonian, we obtain the following equa¬ 


tions of motion 


ddcell 

dt 


dt 


dt 


N 

-y dcell + * E 

(2) 


<b 

1 

pU) 

(3) 

Z5j(t)aceno-^o^+* ^ 

(4) 


where we have included the cavity intensity decay with a 
rate Ki and the spontaneous emission of the atoms with 
a rate 7 . Associated with these decays, are correspond¬ 
ing Langevin noise operators, F’kj, for the cavity decay 
and for the atomic decay [2]. Note, that we have 
neglected dephasing of the atoms, e.g., due to collisions. 
We assume that all the atoms are initially in the ground 
state | 0 ) and that the interaction with the light is a small 
perturbation to the system. We can therefore assume 
that (Tee^ « dll ~ 0. The noise operators describe vac¬ 
uum noise and will never result in either an atomic or 
field excitation. Hence, they will never give rise to clicks 
in the detector (see Fig. [Io](b)) and we can consequently 
ignore them as described in Ref. [2]. Furthermore, we 
treat (Tio(t) as slowly varying in time and formally inte¬ 
grate Eqs. ([^ and ^ to obtain the field operator inside 
the cell 


aceu(i') = - b E / di" / 

^ jzzl do Jo 

X {t")ni (r)<7m ■ (5) 

To find the field at the detector, we need to propagate 
the field through the filter-cavity. The input/output re¬ 
lations for the filter-cavity are 


ddfjiter 

dt 


2 


dfilter “f \/S^2^l/2dcell, 


— \/ ^ 2 / 2 dfilterj 


( 6 ) 

(7) 


with K 2 being the intensity decay rate of the filter-cavity, 
dfiiter describes the field inside the filter-cavity and d de¬ 
scribes the field at the detector. We have again neglected 
any input noise from the cavity decay since it never gives 
a click in our detector and we have also neglected intra¬ 
cavity losses. Formally integrating Eq. (§, and using 
Eq. (Q, gives 

= ( 8 ) 

i=i 

where 

Ojit) = f dt' [ dt" [ 

Jo Jo Jo 


( 9 ) 
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The efficiency is defined as the probability of hav¬ 
ing stored a single excitation in the symmetric Dicke 
state upon detection of a quantum photon. Neglecting 
hi gher order excitations, the atomic state is projected to 
^/p{t)\tf}a{t)) = /(0|a(t)|00.. .0)|0)/ when the quan¬ 

tum photon is detected at time t. Here, |0); is the vacuum 
of the light in the cavity mode and p{t) is the probabil¬ 
ity density of detecting the photon at time t with r] being 
the single photon detection efficiency. Assuming that the 
driving pulse has a duration of tint, the efficiency of the 
write process is 


/o‘°V(OI(V'D|V'a(i))|^df /o”‘|(6»j(f))e|^df 

where we have used Eq. assumed N :$> 1, and have 
defined the ensemble average (.. .)e = jj ■ ■)■ 

To get an expression for Tywritej we need to eval¬ 
uate and {\9j{t)\'^). We now explicitly in¬ 

clude the spatial dependence of the couplings assum¬ 
ing that gj{t) = giy (t) sin{kqZj{t)) and = 

^xy(t) sin{kcZj(t)), where 


— x- (t)-y- (t) 

( 11 ) 

— —y?(t) 

( 12 ) 


with w being the waist of the beams and {xj , yj , zj) is the 
position of the j’th atom. The transverse ccy-dependence 
is assumed to be Gaussian while the z dependence is si¬ 
nusoidal due to the standing wave in the cavity, kg (kc) 
is the wave vector associated with the quantum photon 
(classical field). We have neglected additional geomet¬ 
ric phases in the gaussian couplings since we are always 
considering the product g*j^j- 

The cavity field is a standing wave along the z-direction 
and both modes are assumed to have a node at the cen¬ 
ter of the cell at z = 0. This geometry ensures an ideal 
overlap between the two modes at the position of the mi¬ 
crocell at the center of the cavity. Away from the center 
of the cavity, the overlap of the two modes is degraded, 
which can lead to a detrimental phase difference if the 
length of the microcell along the cavity axis is to long. 
We estimate that for the Cs-cells used in the proof-of- 
principle experiment, the write efficiency is only degraded 
by a factor ~ 0.97 for a cell length of ^ 1 cm |30) . 

To suppress the effect of Doppler broadening of the 
atomic levels, A will be in the GHz range whereas the 
transverse waist of the beam will be on the order 50/rm. 
Consequently, the transverse coupling and the velocity of 
an atom can be considered constant for the integration 
over t'" appearing in 9j{t), which will have a typical time 
scale of 1/ |A — * 7 |. The z-dependence of the coupling, 
however, varies rapidly due to the standing wave in the 
cavity and cannot be assumed to be constant. Writing 
= Zj{0) +vi^h"', where vi^'’ is the z-component 
of the velocity of the fth atom, we can thus perform 
the integration over t'" and adiabatically eliminate the 


optical coherence since we are far detuned. To obtain 
an expression for \{dj{t))\^, we assume that the spatial 
distribution of the atoms is uniform and that the velocity 
distribution of the atoms follows the Maxwell-Boltzmann 
distribution with temperature T. Both distributions are 
assumed to be independent of time. In our analytical 
calculations, we also assume k^ ^ kg = k and that kL^ ^ 
1 such that « 0, but in our numerical simulations 

presented below, we set the difference between kc and kg 
corresponding to the real level structure of Gs where a 
splitting between |0) and |1) is 9.2 GHz. Here, is the 
length of the cell in the beam direction. 

With these assumptions, we obtain 


mt))e 


^ A [(A + ij/2)/rd] 


1 

K 1 K 2 


, (13) 


where we have assumed that e ('^ 1 / 2 )* ~ g (^. 2 / 2)1 ^ 
Furthermore, we have assumed that the cell dimensions 
{x xy X z) are 2Lx2Lx 2L^ and that erf {y/2L/w) « 1 
meaning that we ignore any small portion of the beam, 
which is outside the cell. w[...] is the Faddeeva func¬ 
tion defined as w[z] = (1 — erf(—iz)) and F^ = 

a/2A:bT /mk is the Doppler width of the atomic levels 
at the temperature T where m is the atomic mass and 
fee is the Boltzmann constant. 

While \{dj{t))\'^ does not contain any correlations be¬ 
tween an atom’s position at different times, {\9j{t)\ ) 
does. Eq. (101 thus characterizes the effect of the ran¬ 
dom atomic motion and the motional averaging associ¬ 
ated with it. We evaluate | {9j under similar assump¬ 
tions for the atoms as before. The correlations decay 
in time such that after many collisions with the walls, 
an atom’s position is completely uncorrelated from its 
initial position. We assume that the decay of the cor¬ 


relations is exponential such that, e.g., {gxy i0)gxy (t)) = 
{{gxy {Q))‘^)e~^* + {9xy (0))^(1—e“^*), where the first term 
contains the short time correlations while the second 
term characterizes the long time limit where the correla¬ 
tions are only through the average values. We substan¬ 
tiate this assumption by simulating a box of randomly 
moving, non-interacting atoms and find good agreement 
with a decay rate F = ai^thermai/w where r^thermai is the 
average thermal velocity of the atoms, w is the waist of 
the Gaussian cavity mode and a is a numerical constant 
on the order of unity |30] . Employing this model for 
the atomic correlations and assuming K 2 tint ^ 1 such 
that the effective interaction time (I/K 2 ) is set by the 
linewidth of the filter-cavity, we find that 


-I- \g\^\fl\^ A{m,K2,T) X 

( /Re{vif [(A-I-f7/2)/rd]} 

1^32FdLA ^/2 

^ Im{w[(A-|-^7/2)/rd]} ^ 
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+ ^|z.[(A + z7/2)/r,]|^|^^, (14) 

where we have defined 

4/' _ (2r + Ki+/C2) 

A{k^,K2,T) - + Ki)(2r + /C2 )(ki + /C2) ’ ^ ^ 


and we have neglected all terms oc since these av¬ 

erage to zero rapidly. Using Eqs. (131 and (14), we can 
directly evaluate ? 7 write from Eq. (10). In the limit of 
Ki > (r, K 2 ) and A > Ed » 7 , the expression for ?7write 
reduces to 


^write 


1 


1 I k .2 ( 4 L^ 

2T + K2 \7TW'^ 


1 ) 


1 - 


1 


N, 


pass 


(16) 


where we have assumed L > w. Equation (16) shows 
that 7 write 1 as K 2 /r —>■ 0 , i.e., the write efficiency 
improves with the length of the effective interaction time. 
This is the motional averaging of the atomic interaction 
with the light. Eq. also shows how the efficiency 
improves as the ratio between the beam area and the cell 
area {ttw^/L'^) increases. The last equality in Eq. (16) is 
obtained assuming that T >> K 2 . In this limit Apass ~ 


can be interpreted as the average number of passes 
of an atom through the beam during the decay time of 
the filter-cavity. 

To describe the write efficiency quantitatively, we have 
numerically simulated the experiment with Cs-cells in¬ 
cluding the full level structure of the atoms as described 
in Ref. [30] . The A-scheme level structure can be realized 
with the two ground states |0) = |F = 4,mi? = 4) and 
| 1 ) = |F = 3,tof = 3) in the 6 ^ 51/2 the ground state 
manifold and the excited state |e) = \F' = 4:,mF' = 4) 
in the excited Q^P 3/2 manifold. Note that with this con¬ 
figuration, the quantum and classical field differ both in 
polarisation and frequency and the filtering of the quan¬ 
tum photon is thus expected to be easily obtained us¬ 
ing a combination of both polarisation filtering and the 
filter-cavity. Fig, [da] shows the simulated write efficiency 
as a function of K 2 - It is seen that rj^rite ~ 90% for 
K 2 « 27r • 10 kHz, which translates into a write time of 
~ 150 /iS. Furthermore, we estimate that the number 
of classical photons, which should be filtered from the 
quantum photon is ~ 4.4 • 10® for realistic experimental 
parameters |30|. This level of filtering is expected to be 
easily achieved using frequency filtering. 




FIG. 3. (a) Write efficiency as a function of the linewidth 
of the fifter-cavity. We have simulated a Cs-cell with L = 
150 /rm and w = 55 corresponding to the cells being 
used in the proof-of-principle experiment. We have assumed 
a detuning of A ~ 27r-900 MHz, a pulse length of tint = IO/K 2 
and tvi = 27r • 46 MHz. (b) Optimal readout efficiency as a 
function of the readout time Tread without the filter cavity 
(corresponding to ac 2 —>■ 00 ). The efficiency was simulated for 
the same Cs-cells as the write efficiency and we have assumed 
that Tread = S/Fread where Fjead IS the readout rate, which 
is proportional to the classical drive intensity. The optical 
depth was assumed to be 168 as measured in the experiment. 
The finesse of the filter cavity was varied between 20 and 100 
to get the optimal readout efficiency. We have included the 
full level structure of ^®®Cs in the simulations |30| . 


II. PROOF-OF-PRINCIPLE EXPERIMENT 

To confirm the validity of the model and the results ob¬ 
tained above, we have performed a proof-of-principle ex¬ 
periment, which confirms the most important prediction, 
the presence of a spectrally narrow coherent peak of the 
scattered light arising from motional averaging. While 
several previous experiments [HHH M have demon¬ 
strated long coherence times of room temperature atoms. 


here we wish to show a long coherence time of the emitted 
photons thus demonstrating that the motional averag¬ 
ing technique can be exploited to make coherent photon 
emission. To do this, we compare the theoretical predic¬ 
tions with the experimentally observed power spectral 
density (PSD) of light scattered by the atoms. In this 
proof-of-principle experiment, linearly polarized probe 
light, off-resonant from the atomic transition, interacts 
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(a) 



FIG. 4. (a) Schematic of the proof-of-principle experiment. 
BS: Beam splitter; PBS: Polarizing Beam Splitter; CM: Cav¬ 
ity mirror; PD: Photodiode (b) The D2 transition probed 
in the proof-of-principle experiment, (c) Effective coupling 
scheme for the Faraday interaction [^. The strong probe 
beam (straight arrows) is polarized perpendicular to the ap¬ 
plied field and can thus drive (t+ and (t_ transition. An atom 
scattered between two different mp levels results in a tt po¬ 
larized photon (wiggly lines), orthogonal to the drive. In the 
weak probing limit, the measurement of the Faraday rotation 
angle is thus equivalent to a heterodyne measurement of the 
emitted light in the Raman transition with the probe pulse 
acting as a local oscillator. 


with the atoms resulting in Faraday paramagnetic ro¬ 
tation of the light polarization [5] and the polarization 
state of the light is recorded with balanced polarimetry. 
As explained below, the balanced polarimetry establishes 
a heterodyne measurement of the Raman scattered pho¬ 
tons allowing us to determine their spectrum. 

The experimental setup is shown in Fig. The light 
is red-detuned by 2.8 GHz from the F = 4 i—>■ F' = 5 
D2 transition. For this detuning, the polarisation state 
of the probe light is affected by Cs atoms in both 
F = 4 and F = 3 ground state manifolds. The atomic 
ensemble is contained in a glass-cell with dimensions 
300 /rmx300 /xmx I cm corresponding to an average wall- 
to-wall time of flight of ^ 1.4 /is. The walls of the cell 
are covered with an alkene coating El HZ], resulting 
in longitudinal and transverse spin lifetime in the dark 
Ti « 17 ms and T2 « 10 ms, respectively. The atomic 
density inside the cell is estimated to be ^ 8x 10“^° cm“^ 
[25] . The cell is placed inside a standing wave optical cav¬ 
ity to enhance the light-matter interaction. The cavity 


has finesse F « 17, determined by the output coupler 
(intensity reflection R 2 ~ 80%) and the optical losses in 
the cell; currently the light intensity loss in the cell is 
~ 13% per round trip, limited by the deterioration of the 
anti-reflection coating of the walls during the cell fabrica¬ 
tion. A Pound-Drever-Hall technique is used to lock the 
cavity on resonance. The cavity mode has « 55 /rm waist 
radius, which is a compromise between the requirement 
for strong coupling of light to the atomic ensemble and 
the requirement for low propagation losses through the 
cell. A small portion of the beam at the cavity output 
is used in a feedback loop to compensate for the probe- 
intensity drift and maintain the same photon shot noise 
during the time of measurement. 

A DC bias magnetic field perpendicular to the probe 
direction sets the Larmor frequency of the atoms. Due to 
technical limitations related with the phase noise of the 
laser and the cell birefringence, the polarisation of the 
probe was at an angle of approximately 40-45 degrees 
with respect to the axis of the magnetic field. When 
the probe light is far detuned, the Faraday rotation is, 
however, independent of this angle [5|. For simplicity, 
we therefore describe the dynamics using the level struc¬ 
ture in Fig. ^c), which assumes that the driving field is 
CT+ -I- CT_ polarized, perpendicular to the direction of the 
magnetic field tt. In the far detuned limit, the Faraday 
rotation is due to Raman transitions between magnetic 
states with magnetic quantum numbers uif differing by 
±1. In these Raman transitions, a 7r-polarized photon is 
emitted as shown in Fig. ^c). In the balanced polarime¬ 
try, the driving field and the scattered tt component of 
the light are mixed on a polarising beam splitter and 
the difference intensity is recorded. This corresponds to 
the driving field acting as local oscillator for a hetero¬ 
dyne measurement of the emitted 7r-polarised light. The 
recorded spin noise is thus a measurement of the pho¬ 
tons emitted from the atoms through Raman scattering 
between the Zeeman sublevels of the Cs hyperfine man¬ 
ifolds and is therefore exactly the quantity we are inter¬ 
ested in for probing the coherence of the emitted photons 
and verifying the predictions of the model. The exper¬ 
iment is performed in the continuous regime with con¬ 
stant laser intensity. By comparing the emitted light to 
the shot noise level, we can extract the Raman scattering 
rate. For a pulse duration that can lead to an efficient 
write step, e.g. t = 106 /rs, corresponding to «;2 = 27r • 15 
kHz in Fig. [^a), we find that ^8 Raman photons are 
scattered in the upper-sideband mode EOj. Due to the 
linearity of the process, the spectrum is expected to be 
the same at the single photon level. 

The measurement is performed on atoms in approxi¬ 
mately their thermal state i.e., the atoms are randomly 
distributed in the 16 m-sublevels of the F=3 and F=4 
hyperfine manifold. There is a small deviation from the 
thermal state due to weak optical pumping from the 
probe, and all measurements are recorded in the resulting 
steady state. In this case, there is no macroscopic orien¬ 
tation and the probe induced back-action noise is negligi- 
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ble. The polarimetry noise is the sum of the photon shot 
noise and the Raman scattered photon noise. The photon 
shot noise has a white power spectrum, whereas the spec¬ 
trum of the recorded Raman noise is centered around the 
Larmor frequency due to the energy difference between 
magnetic sub-levels. We perform spin noise measure¬ 
ments for two different Larmor frequencies, ~ 0.8 MHz 
and ~ 2.6 MHz, at the same probe power. By subtract¬ 
ing the two power spectra, the photon shot noise and the 
electronic noise contribution to the recorded spectra can 
be removed and the Raman noise is acquired. 

The measured PSD and the simulations of it are shown 
in Fig. The measured Raman noise reflects two differ¬ 
ent correlation decay timescales: a fast decay timescale 
^ 1 ^s associated with the transient time of flight through 
the probe beam, and a relatively slow decay ^ 100 /xs, 
due to the spin decoherence (probe induced spin relax¬ 
ation). Since the spectrum shown in Fig. [^is measured 
for the scattered light, it is seen that the long spin coher¬ 
ence translates into a long coherence time of photons at 
the single photon level consistent with the theory. The 
PSD is recorded with a higher frequency resolution than 
shown in the figure but we bin the data with a frequency 
resolution (A/ « 61 kHz) chosen so that the spin noise of 
atoms in the two hyperfine manifolds associated with the 
slow correlation decay timescale is contained in a single 
frequency bin. By doing this, complications arising from 
the nonlinear Zeeman splitting and the difference in the 
gyromagnetic ratio between the different hyperfine man¬ 
ifolds can be ignored. The simulations were carried out 
as described in Ref. [30] and have been rescaled to fit the 
measured SPD at f=823.7 kHz. We have carried out sim¬ 
ulations both where the atomic collisions with the wall 
coatings happen instantaneously (ttrap = 0 ), i.e. so that 
the trapping time is negligible compared to the transient 
time, and with a trapping time of ttrap — 0.1 ^s. Fig. 
shows an excellent agreement between the experiment 
and the model with zero trapping time. From Fig. 
we estimate that any trapping time in the experiment 
is ^ 0.1 ^s and can thus be ignored. The narrow peak 
in the scattered light is due to the fact that atoms re¬ 
peatedly come back into the beam with the same spin 
phase, whereas the broad background is due to single 
transients through the beam. The narrow peak in the 
data thus demonstrates the long coherence time of the 
forward scattered light due to motional averaging. Con¬ 
sidering the random motion of the atoms, the only coher¬ 
ence that can survive for that long is linked to the sym¬ 
metric Dicke state as described in the theory. In essence, 
the idea of the motional averaging is to use a spectrally 
narrow filter cavity to select only the photons emitted in 
the narrow coherent peak. Since this narrow peak cor¬ 
responds to a long interaction time, this means that all 
atoms participate equally in the resulting spin wave. Fur¬ 
thermore, since the narrow peak is much higher than the 
broad background, the loss in efficiency from the spectral 
filtering is limited. The excellent agreement between the 
simulation and the experiment thus confirms the applica¬ 


bility of the motional averaging as well as the theoretical 
model we use. 


III. READOUT 

We now consider the readout process. Assuming a sin¬ 
gle excitation has been stored in the symmetric mode in 
the ensemble, a classical drive (D) is applied to read out 
the excitation as a cavity photon. The relevant Hamilto¬ 
nian is obtained by interchanging and in Eq. Q . 
From Heisenberg’s equations of motion, we obtain 


ddcell 

dt 




N 


- 0 ) 

'Oe 


(17) 


i=i 


= - (^ - ^Oe + i9jit)aceii + (18) 


dt 


da. 


U) 

01 


dt 


= i- 


mt) 


‘O') 

^Oe ’ 


(19) 


where we have assumed that dee —^oo ~ and that the 
dynamics of a[{)^ are governed by the classical drive (fl). 
Furthermore, we have neglected the noise operators as¬ 
sociated with spontaneous emission and cavity decay, as 
in the write process. We can formally integrate Eq. (18 1 , 


assuming the xy-dependence of the couplings to be con¬ 
stant for the integration while the z-dependent parts are 
of the form sin(fc( 0 j(O) -I- vi^\o)t)), as in the write pro¬ 
cess. The integration gives a set of coupled equations 


where 


^ (j) 

= ^3 (i)aceU + Cj (t)dlfj ^, 


U) 

01 


( 20 ) 

( 21 ) 


A(t) = 


1 ^ , 


i=i 


/,( 7 ,A,A:) ( 22 ) 




xy ' 


= Ye /i(7,A,A:), 

g2ikzj{t) _ ^ 


/i(7, A,A:) = 


^ 7/2 - i{A - kvi^\t)) 

1 _ g-2ifczj(t) 

7/2 - f(A -I- kvi^\t)) ^ 


(23) 

(24) 


(25) 


and we have assumed that kc ~ k, 
can be expressed as a matrix system of the form 


q ~ k. The equations 


dx(t) 

dt 


= M(t)x(t), 


(26) 
















where x(t) = ^Oceii, • ■ ■ > ^(0 is the 

coupling matrix between the atoms and the cavity field. 
We now split the couplings into (large) average time- 
independent parts, {(A}e, (^}e} = {A,B,C}, and 

(small) time-dependent parts, {SA{t),SBj{t),SCj{t)}. 
The coupling matrix can then be expressed as M(t) = 
Mo -I- (5M(t) where Mq contains the average time- 
independent couplings while (5M(t) contains the time 
dependent fluctuations. Assuming the readout pulse is 
long, the atoms will have had the same average interac¬ 
tion with the light meaning that M(t) « Mq. Treating 
(5M(t) as a small perturbation, we can then obtain a per¬ 
turbative expansion of Sceii- Assuming that the initial 
state of the atoms before readout is the symmetric Dicke 
state, we find that, to second order in (5M(f), the cavity 

field can be expressed as ficeii ~ ®c2i + ®ceii- Here we have 
omitted the first order term, which depends on the fluc¬ 
tuations in dA{t), since we find that its contribution is 
suppressed by a factor of at least dF/N compared to the 
other terms. Here dlTca.-v oc Ng ^is the optical depth 
on the |0) e-)- |e) transition per cavity roundtrip time and 
F = 27r/(TcavKi) is the finesse of the cell-cavity. 

As in the write process, the field from the cell-cavity 
is sent through the filter-cavity in order to both filter 
the classical drive photons from the single photon and to 
filter out incoherent photons as we will describe below. 
Using Eqs. @- 0 , we find that the readout efficiency is 

4 Jo Jq 

x(aLn(i')Sceu(t")), (27) 


where Tj-ead is the duration of the readout stage. To lowest 
order we find 


^cell {i) 


( 0 ) ^ VNB _ ^-Fvt 


“cell “ 


Vv 


(1 - , (28) 


where V = {C — A)^ + ANB'^. Inserting Eq. (281 in 
Eq. (271 and taking the limit of H —)■ 0 and Tread j K 2 00 , 


gives a zeroth-order readout efficiency of 


f?read,0 ~ ^ , i ’ (29) 

djr “T 

Eq. ( [ 2 ^ is equivalent to the result for cold atomic en¬ 
sembles [2] and represents the long time limit of perfect 
motional averaging where the efficiency improves with 
optical depth and finesse of the system. 

The coherence time of real atoms is, however, limited 
and a fast readout is therefore desirable. Erom Eq. ( [^ , 
we identify the readout rate Tread = \{A+C-\-^/V). This 
readout rate Tread increases with increasing strength of 
the readout pulse, so that for strong driving correspond¬ 
ing to a fast readout, it is necessary to consider higher 
order terms in the perturbative expansion of aceii(t)- To 
second order, we find that ?7read ??read,o + ??read ,2 where 
the second order term (77read,2) mainly describes the loss 


of the excitation due to spontaneous emission. Con¬ 
sequently, the magnitude of ? 7 read ,2 increases with the 
driving strength while it’s sign is negative. ? 7 read ,2 con¬ 
tains correlations between an atom’s position at differ¬ 
ent times, which we can treat in a similar manner as in 
the write process, i.e. as exponentially decaying in time. 
By simulating the readout process with the Cs-cells in a 
similar fashion as for the write process, we can quantita¬ 
tively describe the readout efficiency to second order m- 
Eig. 1 ^ shows the readout efficiency to second order as a 
function of the readout time. We have assumed an opti¬ 
cal depth of 168 and varied the finesse of the cell-cavity 
between 20-100 to get the maximum readout efficiency. 
The readout time is set to T^ead = S/Tread ensuring that 
a negligible population is left in the system at the end of 
the read out stage. (Note that in these simulation we do 
not include the filter cavity considered for the write stage. 
Eormally this corresponds to taking the K 2 —t oo). The 
full level structure of ^^^Cs is included in the simulations 
and the optimization in the finesse is due to the extra lev¬ 
els in Cs-atoms, which introduces additional couplings. 
In general, high (low) finesse is optimal for short (long) 
readout times. A small cavity detuning was also included 
in the optimization in order to compensate for the shifts 
caused by the couplings to the extra levels [30] . For a fi¬ 
nesse of ^50 and a readout time of t « 200 gs, a readout 
efficiency of reread ~ 90% is obtained. 


IV. ERRORS 

Above we have focussed on the efficiency of the pro¬ 
tocol. We will now consider the errors, which limit the 
performance of the system as a single photon source with 
memory. We find that the dominant errors are multiple 
excitations during the write process and the possibility of 
reading out atoms, which have been incoherently moved 
to state |1) by either inefficient optical pumping or wall 
collisions. 

Multiple excitations in the write process would also 
create multiple quantum photons, which could in prin¬ 
ciple be discriminated from the situation with a single 
quantum photon if perfect single photon detection is pos¬ 
sible. In a realistic setup there will, however, always be 
some finite detection probability r]^ and the probability 
of creating two excitations would introduce an error of 
~ 2(1 — ??d)Pe to lowest order where Pe oc /g‘"*(|0j(t)|^) is 
the excitation probability. This error can be made arbi¬ 
trarily small by simply decreasing pe, i-e. decreasing the 
strength of the classical drive. This will, however, also 
decrease the rate of the operation, which scales as 1/pe- 

Atoms can also be in the readout state |1) either by 
inefficient optical pumping or by wall collisions. These 
atoms will mainly produce incoherent photons different 
in both frequency and temporal shape from the coherent 
single photons originating from the symmetric excitation. 
The incoherent photons will have a much broader tem¬ 
poral and frequency profile than the coherent photons. 
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FIG. 5. The probability to read out incoherent photons (pi) 
normalized by the fraction of atoms (e) that have been inco¬ 
herently transferred to the readout state (11)). Pi /e essentially 
only depends on Tread^read and the linewidth of the filter- 
cavity (k 2 ) for the parameters that we are considering, which 
are e <C 1, an optical depth of 168, and a finesse of the cell- 
cavity in the range 20-100. Furthermore, we have assumed 
that Tread^read = 3 which eusures a temporal filtering of the 
incoherent photons while keeping a high readout efficiency of 
the coherent photons. The plot was obtained by numerically 
simulating the Cs-cells used in the proof-of-principle experi¬ 
ment including the full level structure of the Cs-atoms. 


We can thus to some extent filter them from the coher¬ 
ent photons by sending the light through a filter-cavity, 
which makes a spectral filtering, as well as having a not 
too long readout time Tread, which makes a temporal fil¬ 
ter. In addition to the incoherent photons, atoms inco¬ 
herently prepared in the wrong state can also produce 
coherent photons because the incoherent atoms have an 
overlap with the symmetric modes. If a fraction e of the 
atoms are transferred to the state |1), the probability 
to read out a coherent single photon from these atoms 
is epread- The probability pi to read out an incoherent 
photon can be found to lowest order by assuming that 
an excitation is stored in any asymmetric mode instead 
of the symmetric Dicke mode in the perturbative expan¬ 
sion of CLceii described above. Doing the perturbative ex¬ 
pansion, we then get a contribution to Oceii from these 
incoherent excitations to the first order term oi. From 
this, we can find the number of incoherent photons in the 
retrieval. We have evaluated pi by numerical simulating 
the Cs-cells as for the readout n. Fig. I shows Pi/e as 
a function of the linewidth (^ 2 ) of the filter-cavity. We 
have assumed that T^ad = S/Tgead as in Fig. 3b Note 
that this choice of readout time ensures a high readout 
efficiency while still making a temporal filtering of the in¬ 
coherent photons since these have a smaller readout rate 
and hence predominantly arrive later. It is seen that 
Pi « e for K 2 ~ Stt • 80 kHz. With a linewidth of the filter 


cavity more narrow than this, the error will thus be dom¬ 
inated by the coherent photons which are emitted with a 
probability erj. Imposing this linewidth of the filter cav¬ 
ity for the numerical example for the readout efficiency 
given at the end of Sec. Ill with a readout time of f = 200 
ps would make it drop from « 90% to « 86%. Hence we 
loose only a little on the readout efhciency by filtering 
out the incoherent photons. Experimentally, it will be 
simpler to use the same filter cavity for the retrieval as 
for the write process, and hence it may be desirable to 
use a more narrow filter cavity to have an efficient write 
process (see Fig. 3b). In this case one can use a longer 
read out time r^ead to suppress loss from the filter cavity. 
After filtering out the incoherent photons, the remaining 
error is thus only due the coherent photons from atoms 
being incoherently prepared in the wrong state, and the 
error is equal to the probability e that each atom is in 
the wrong state. 


V. CONCLUSION AND DISCUSSION 

In conclusion, we have developed a theory for mo¬ 
tional averaging for discrete variable systems and pro¬ 
posed an efficient and scalable single photon source based 
on atomic ensembles at room temperature. We have con¬ 
sidered a specific setup where the atomic ensemble is kept 
in a small cell inside a cavity and shown how both read 
and write efficiencies above 90% can be achieved for a 
real experimental system based on Cs-atoms. The write 
and read processes have a timescale of 100—200 /iS, which 
is considerably shorter than the demonstrated quantum 
memory time of 10 ms [28]. To verify the essential effect 
described by the theory, we have performed a proof-of- 
principle experiment with room temperature Cs atoms 
contained in a microcell with spin preserving coating de¬ 
posited on the walls. The measurement of the scattered 
light reveals long coherence time at the single photon 
level, resulting in a narrow peak, which is in excellent 
agreement with the theoretical model being used. This 
thus confirms the essential feature of the theory. 

The room temperature cells considered here provides 
a promising building block for future quantum networks 
because of their scalability compared to cold atomic en¬ 
sembles. As a particular application, we have considered 
a basic step of DLCZ quantum repeater with a single 
entanglement swap and a distance of 80 km assuming 
a dark count rate of 1 Hz and single photon detection 
efficiency of 95% jH Hj. Including various experimental 
imperfections such as intra-cavity losses and inefficient 
in/out coupling of the cavities, we estimate that a pair 
with fidelity ~ 80% can be distributed with a rate of 
~ 0.2 Hz [30]. In this estimate, we have neglected effects 
from limited memory time and have assumed that a frac¬ 
tion of 0.5% of the atoms have been incoherently trans¬ 
ferred to the state |I). Note, however, that the rate of 
entanglement distribution can be greatly enhanced using 
spatially multiplexing schemes, which are possible due to 
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the scalable nature of the room temperature cells. A par¬ 
ticularly attractive feature of such multiplexing is that it 
also decreases the necessary memory time cni, and thus 
relaxes one of the most challenging requirements for long 
distance communication based on atomic ensembles. The 
microcells introduced here may thus serve as an essen¬ 
tial building block for future photonic networks. On the 
other hand, for more near term applications the scalable 
nature of the setup will also be highly interesting for ap¬ 
plications requiring multiple single photon inputs such as 


for instance photonic quantum simulators [IMI]- 
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VI. SUPPLEMENTAL MATERIAL: SCALABLE PHOTONIC NETWORK ARCHITECTURE BASED 
ON MOTIONAL AVERAGING IN ROOM TEMPERATURE GAS 


In this supplemental material to the article “Scalable photonic network architecture based on motional averaging 
in room temperature gas”, we describe the details of the numerical simulations referenced in the article. We also 
estimate the number of classical photons, which should be filtered from the quantum photon and present the exact 
definition of the optical depth mentioned in the article. Finally, we estimate the errors limiting the performance of 
the proposed setup in a DLCZ-like quantum repeater. 


VH. NUMERIGAL SIMULATION - WRITE 


To justify our assumption of an exponential decay of the correlations appearing in (|dj(t)|)^ and to qualitatively 
characterize the write efficiency, we perform a numerical simulation of a gas of non-interacting atoms in a cell. We 
have based the simulation on the microcells filled with Cs-atoms, which were used in the proof-of-principle experiment. 
These cells have dimensions of 300 /rm x 300 /rm x 1 cm. The cells have been placed inside a cavity with a linewidth 
of Ki « 27r • 46 MHz and both the field from the quantum photon and the classical drive are assumed to have 
approximately a Gaussian shape with a waist of 55 ^m. The small beam waist ensures that we can neglect cavity 
losses from the walls of the cell. An approximate A-atom can be realized in the hyperfine states of Cs with state 
|0) = lA = 4, tuf = 4) and state |1) = |F = 3,tof = 3) in the 6^5'i/2 ground state manifold. The Doppler width of 
the atomic levels is T^ 27r • 225 MHz at a temperature of T = 293 K and we assume a detuning of A ^ dT^ from 
the excited level such that the effect of Doppler broadening is negligible (see below). 

Starting from Eq. (9) in the article and performing the integral over t" as described in the text below, we can 
express {\dj{t)\)‘^ as 


{m)f)e = Y^ dt[ 


dt'l 


dt' I dt^^ € K2l2{t ^ 2 )^ ^i/2(i2 ^ 2 ) 


X {XY* {t'DXY, (t")). 


where we have defined 


zAt) 


XY,{t) 


gi(Afc)zj(t) _ ^-i(k^+kq)zj(t) 


- 7/2 -I- i(A -I- kcvi^\t)) 
rtge ^ . 


g-i(Afc)zj(t) _ gi(fcc+fc,)zj(t) 

- 7/2 -I- i{A - kcvi^\t)) 


(30) 


(31) 

(32) 


Note that we have not made the assumption oi k^ ^ kg Ri k, as in our analytical calculations, since we have the 
27r • 9.2 GHz splitting between the ground states, which corresponds to kg — kc = A^ « 193 m“^. With a cell length 
of 2Lz = 1 cm the assumption 2 AfcLz -C 1 is close to being violated and we shall therefore not make this assumption. 
The extra terms oc will approximately result in a factor of 


_ (cos (AfcZj-))^ 
(cos(AfcZj)^) ’ 


(33) 


which should be multiplied with the analytical expression for the write efficiency, which was obtained assuming A^ = 0 . 
Fig. [ 6 ] shows how CAfc depends on the length of the cell assuming that the atoms are equally distributed in the entire 
cell. It is seen from Fig. |^that as long as the length of the cell is 2Lz < 1 cm then CAfc ^ 0.97 for Aj, « 193 m“^ 
and hence the frequency difference between the quantum and classical fields does not significantly degrade the write 
efficiency. In all our numerical simulations, we, however, keep the terms oc for completeness. 
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FIG. 6. Limitation in efficiency, CAfe, from the difference in wavenumbers as a function of the cell length 2Lz for Ak ~ 193 m”'^ 
corresponding to the 27r • 9.2 GHz splitting between the hyperhne ground states of ^^^Cs. The atoms are assumed to be evenly 
distributed in the cell, ca^ > 0.97 for 2Lz < 1 cm. 


The correlations appearing in {XY*{t'()XYj{t2)Z*{t'l)Zj{t2))e (see Eq. (30)) depend on the time difference \t'( — t2\ 
and we therefore introduce the shorthand notation {XY*{t'()XYj{t2)Z* {t^JZj{t2))e = {XY, Z)e{t'{ — t',£j. We change 
to the variables u = t'{ +and s = t'( — and by changing the order of integration, we can perform the integrals over 
t,t'i,t2 and u. To obtain the write efficiency rywritej we need to perform an additional integration over t (see Eq. (10) 
in the article). We are therefore left with 


{\0medt = 


h{tint,Ki, K2, s){XY, Z)eis)ds, 


(34) 


where ki, K2, s) is a function of s obtained by performing the integrals over t,t'i,t'2 and u. We can evaluate 

the integral over s numerically by simulating the correlations {XY,Z)f,{s). Since the atoms do not interact with 
each other, we independently simulate the motion of fV = 5000 atoms through the cell and evaluate the correlations 
of atoms at points separated in time by s. Finally, we average over many realizations. At room temperature we 
assume that the motion of the atoms is classical. Furthermore, we assume that the atoms are evenly distributed in 
the cell and that their velocity distribution is described by the Maxwell Boltzmann distribution at a temperature 
T = 293 K. We assume that the atoms are re-thermalized completely after every collision with the walls of the cell 
but qualitatively similar results are obtained for a ballistic model without thermalization. For the ballistic model, 
the Zj parts of the couplings in principle do not average down. This could lead to effects not averaged away by using 
narrow filter cavities. However, we are far detuned compared to the Doppler width of the atoms and the cavity fields 
are standing waves, which can be viewed as the superposition of two counter propagating waves. As a result, the 
effect of the velocity fluctuations of the atoms cancel and the fluctuations in the Zj terms are greatly suppressed. 
This is in contrast to what happens in ensemble based schemes with a laser coming from one side where Doppler 
effects do not go away by working far off resonance [5]. Consequently, we expect similar results for the ballistic 
model as for the model with complete thermalization. We have explicitly verified this by repeating our simulations 
with the ballistic model, which lead to more or less identical results as for the model with thermalization (the two 
simulations cannot be distinguished due to the noise from the random nature of the simulations). The result of a 
simulation with thermalization is seen in Fig. |8a[ which shows how the correlations decay as a function of s such that 
for s —>■ oo, we have (XY,Z)e{s) —)■ \{XY)e\ |(^)e| ■ This enables us to introduce a maximal cutoff, Smax, in the 
numerical integral appearing in Eq. (34), above which, the correlations have effectively vanished. As a result, we can 


semianalytically evaluate rewrite for an arbitrary pulse length tjut without additional numerical difficulty. Note that 
Fig. 1^ also shows that the exponential model of the decay of the correlations assumed in our analytical calculations 
is a good approximation. 

Based on our simulations of {XY, Z)e[s), we have also estimated the power spectral density, PSD in the proof of 
principle experiment. The PSD was measured (see Fig. 2 in the article and Fig. below) by measuring the Faraday 
rotation of light transmitted through the cell. In these measurements the light is emitted by the scattering between 
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different mp states in an applied magnetic field. As a result, the signal is modulated at the Larmor precession 
frequency. We therefore consider the modulated power spectral density 

1 piint ptint 

PSD(/)oc— / dt 

^int Jo Jo 

1 pJ'int p^int 

OC^/ dt 
^int Jo Jo 

1 p^int ^^int 

+ 7 ^/ dt dt'i(XY,Z},it-t')-\{XY),\^\{Z},\^)e^^^fd-tJ 

*int Jo Jo 

1 p^int ptint 

^df,o + ^ dt dt\{XY,Z),{t-t')-\{XY),f\{Z),f)e^^^dd-t') 

^int Jo Jo 

where / is the frequency and Sf^Q is the Kronecker delta function. Note that we have assumed that the coherent 
light-atom interaction gives a contribution of ^ /*'“* dt J*'"* dt' \{XY)e'^ \{Z)e'^ g2*-n-/(t-t) _ pgp), jn the 

int 

proof-of-principle experiment, the atoms are subject to a magnetic field, which makes the atomic spins process around 
the mean spin direction with a Larmor frequency of 823.8 kHz. The Raman scattered part of the PSD is thus centered 
around this frequency. However, the measured PSD also contains both the shot noise of the light and electronic noise 
from the measurement equipment. Since we are only interested in the signal from the atomic interaction, this noise 
is isolated by performing a second measurement at a higher Larmor frequency (2594 kHz) and subtract it from the 
first measurement. The higher Larmor frequency is chosen such that the two atomic signals are well separated in 
frequency. In the simulated PSD, we have fitted a Lorentzian to the broad feature centered at 823.8 kHz and shifted 
it to be centered at 2594 kHz. This has then been subtracted from the simulated data to include the substraction of 
the two signals in the experiment. 

To validate that the proof-of-principle experiment is probing the theory in the right limit of single photon Raman 
scattering, we estimate the number of Raman photons scattered over the relevant pulse length. This number can be 
found from the ratio of the height of the central peak in the PSD to the shot noise of light. In the experiment, the 
balanced polarimetry measures the amplitude quadrature of the field in the polarisation axis perpendicular to the 
drive field dsc- Neglecting proportionality factors irrelevant for this calculation, the PSD at a discrete frequency / is 
given by: PSD(/) ex AI(/)A(/), where: 


(35) 

(36) 

(37) 


-^(Z) oc / [dscit) + a|j,(t)] e^'^f^dt = dsdf) + [asc(-/)]^ • 


(38) 


Here dsdf) = /q™ dsdt)^'^™^dt and t^s is the measurement time. It can then be shown that the PSD takes the form: 


PSD(/) oc I -k [dsc(-/)]^ [asc(-/)] + [asc(/)]^ [dsdf)] ■ 


(39) 


The unity term in the left-hand part of the above relationship is the photon shot noise contribution, present even when 
there is no scattered field; the second and third terms represent the Raman-scattered number of photons in the lower 
(Stokes line) and upper (anti-Stokes line) sideband respectively. In Fig. the PSD recorded with a pulse duration of 
tms ~ 16 /iS is plotted. The peak height at the Larmor frequency due to the Raman scattering has equal contributions 
from the Stokes and the anti-Stokes lines. The number of coherently scattered photons during the measurement pulse 
is the ratio of the peak height excluding the photon shot noise and the incoherent photons (denoted with a dotted 
lines in Fig. to the photon shot noise level, which corresponds to a power of one photon per unit bandwidth as 
shown in Eq. ( |39[ ). From Fig. we find that for the proof-of-principle experiment approximately 1.2 photons are 
coherently scattered in the Stokes line during the 16 /rs pulse. This corresponds to approximately 8 photons over a 
pulse duration of 106 /iS that can lead to an efficient write step according to our proposal. Although this is slightly 
higher than the single photon required for the protocol, we expect the spectrum to be independent of intensity at 
these light levels. 

In the above calculation it is important to estimate the photon shot noise level. For this, the electronic noise, 
as measured with no light at the detector, is removed from the recorded PSD and the spectrum is corrected for 
the detector frequency response. The photon shot noise then corresponds to the level where the PSD levels off at 
frequencies a few linewidths of the broad Lorentzian away from the Larmor frequency; at these frequencies the power 
from the Raman scattered photons is negligible and the PSD is determined only by the photon shot noise. The 
photon shot noise is also verified by performing balanced polarimeter detection of light that has not interacted with 
the atomic ensemble and with power equal to the one as the drive light in the experiment. We note that the detector 
response was characterized by such a measurement. 
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FIG. 7. Experimental data of the power spectral density (PSD) of the Raman scattered light measured in the proof-of-principle 
experiment. These data were obtained for a pulse length of ~ 16 ns. The electronic noise has been removed from the data 
and the shot noise has been normalized to unity. The number of coherent Raman-scattered photons, equally distributed in the 
Stokes and anti-Stokes lines, can be estimated by the comparing the height of the central peak (dotted line) to the shot noise 
of light. For the proof of principle experiment, the number of coherently scattered photons in the antiStokes sideband over the 
16 /rs of the pulse is estimated to be ~ 1.2 photons. 


VIII. NUMBER OF PHOTONS 


As mentioned in the article, the purpose of the filter cavity is both to increase the averaging time and to filter the 
quantum photon from the classical photons. We will now estimate the number of classical photons, which needs to 
be filtered from the single quantum photon. In order to do this, we need to characterize the ensemble, which we do 
by introducing the optical depth d. 

To obtain an expression for the optical depth, we assume that we are working with the previously mentioned 
Cs-cells. The relevant level structure is shown in Fig. | 8 b[ In the write process, the classical drive is applied on the 
transitions between the ground state |0) = |A = 4, wf = 4) and the excited states |ei) = \F' = A,mF' = 4) and 
|e 2 ) = \F' = 5,mF' = 4) characterized by gi and ^ 2 , respectively. The quantum photon is created on the transition 
|ei) —>■ |1) = jA = 3,mF = 3) characterized by g. Note that with this field configuration, the cell-cavity in principle 
also mediate the transition jA = 4, tuf = 4) —>■ IF" = 4, ttif =3) in the write setup but this transition is suppressed 
by the 27r • 9.2 GHz splitting between the ground states, which makes the corresponding photon non-resonant with 
the subsequent filter-cavity. This transition will, therefore, never give a click in the detector. Since the interaction is 
only a perturbation to the system, we can therefore neglect this transition in our numerical simulations. All atoms in 
the ensemble are initially pumped to the ground state |0) and, in order to characterize the optical depth, we assume 
that the transition characterized by g is non-driven and ignore any cavity (Purcell) enhancement of the corresponding 
decay. The cavity field thus only couples |0) —>■ |ei) and |0) —>■ |e 2 ) with coupling constants gi and 32 , respectively 
(see Fig. 8 b). 

The equations of motion for the cavity field, Ocav, and the relevant atomic operators in a suitable rotating frame 
are 


dcav = -(K/2)dcav + * + 92\i)K.,Q 

i=i 

^eiO = -{71/2 - + i9iHt)dc^^ 

= -( 72/2 - iA2)d-^elo + i 92 \t)ac^v, 


(40) 

(41) 

(42) 


where = |e/)j(0| {I = 1,2) and we have assumed that ai{li — CTgg^ « —1. For simplicity, we have assumed 
the couplings (g) to be real. Ai = wi — Wcav (^2 = <^2 — ^cav) is the detuning of |ei) (|e 2 )), while 71 ( 72 ) is the 
corresponding decay rate. Here wi ( 0 ^ 2 ) is the frequency associated with the atomic level and Wcav is the frequency of 
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F'=5 

F'=^ 

F'=3 


F=4 

F=3 

(b) 


FIG. 8 . (a) Simulation of the correlations {XY, Z}e{s). The correlations are normalized to be unity for s —^ oo where there 
are no correlations and {XY, Z)s{s) —>■ \{XY)e\^ \ {Z)e\^. The data from the simulation have been fitted with an exponential 
model validating our assumption of an exponential decay of the correlations. The fit gives a decay rate of T = 2n ■ 0.75 MHz 
corresponding to T ~ 1.3uthermai/rc. (b) Sketch of the &^Sx /2 and 6 ^P 3/2 hyperfine levels in ^®^Cs. A A-atom is realized with 
|0) = \F = 4, m_F = 4), |1) = IT = 3,mF ~ 3) as ground states in 6 ^S'i /2 and |ei) = |T' = i, 171 ^ 1 ^ 4 ) as the excited level 
in 6 ^P 3 / 2 - To characterize the optical depth, we consider the transitions |0) —>■ |ei) and |0) —>■ |e 2 ) = |T^ = 5,mFi = 4) 
characterized by gi and § 2 , respectively. 


the cavity field. Formally integrating Eqs. (|41[)-(|42[), assuming that = (Xg^Q = 0 at time t = 0, and inserting the 




resulting expression for o-g'^Q and cTg'^g into Eq. (|40[) gives 


U) ■ 


N 


dcav — (^/2)hcav ^ ^ 

i=i 




o-(7l-*^l)(*-t') -U) 


5l (i)«cav(i')di' 


^0 


(43) 


where we have explicitly written the time dependence of Ocav inside the integrals. We evaluate the integrals in Eq. (43) 
assuming that we can treat aca,v{t') as a constant in time and move it outside the integrals. Furthermore, we assume 
that g\^\t') = sm.(k{zj{Q) + vi^\o)t')) similar to the procedure described in the article. Note that k is the 

wavenumber associated with the cavity field while Zj(0) is the z-part of the position (velocity) of the j’th atom. 

After evaluating the integrals, we obtain 


N 


= -(K/ 2 )ag.. + ^ 5 : 


bqlyW + gi^lyit) Zj{A2,-f2,k) 




(44) 


where we have adiabatically eliminated the optical coherence and have rewritten Zj{t) defined in Eq. (31) to 

g2ikzj{t) _ 2 ^ 2 — 


ZjiA,j,k) = 


7/2 + i{kvj{0) — A) 7/2 — i{kvj{0) + A) ’ 


(45) 


such that 7 , A, and k become variable parameters. We now perform an ensemble average of Eq. (44) assuming that the 
atoms are evenly distributed in the cell and that their velocity distribution follows a Maxwell Boltzmann distribution, 
as previously considered. Furthermore, we assume that the cry-dependence of the couplings are Gaussians similar to 
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Eqs. (11)-(12) in the article and that we are detuned far from the Doppler width of the atoms. This results in 

dcell.^ 


— (/t/2)dcell 


where the imaginary part is contained in 


I 51 IS 1 


+ 


I 52 IS 2 


72/4 + A 2 72/4 + A 2 


8 L2 


(46) 


The second term in Eq. (461 is identified as the single pass optical 


depth, d, divided by the cavity round trip time, r, where exp(—d) is the attenuation of the light field after passing 
through the ensemble. Since d depends on, e.g., the detuning, it is not a direct characterisation of the ensemble. 
Instead, in analogy with Ref. [^, we characterise the ensemble by d, the hypothetical optical depth, which would 
be obtained for resonant fields in the absence of Doppler broadening and hyperfine interaction, i.e., Eq. (46) with 
Ai = A 2 = 0. Furthermore, we assume that 71 = 72 = 7 such that the optical depth is 


d = 


Nt 


51 1 




Otx 


(47) 


where we have defined the factor axy = Note that Eq. (47) can be rewritten to the following well known formula 
for the optical depth [3] 


d = Qtt 


^ ( 'li’O + "I2.0 




(48) 


where 2L is the transverse size of the cell, A = X/2tt is the rescaled wavelength of the light, and 7s_o is the spontaneous 
decay rate of level back to |0) (s=l,2). The optical depth can also be related to the Faraday rotation angle, 9p, 
which is typically measured in experiments and used to estimate the number of atoms, N in the ensemble [SI [B] . For 
^^^Cs, the relation between Op and A is |B] 


N = 


32TrL^epA2 


ai(A2)7A2 


where ai(A 2 ) is the vector polarisability given by 


.AN W 35 

1-A,,5'/A 2 


21 


1 — A 4 / 5 //A 2 


+ 176 


(49) 


(50) 


with Ax'b' denoting the hyperfine splitting between level F' = x and F' = 5. Combining Eq. (48) and Eq. (49) gives 
the following relation between d and 9p 


d = 


12 A 20 _f 7i^o + 72,0 
ai(A 2 ) 72 


(51) 


For the cells used in the proof-of-principle experiment, the Faraday rotation angle has been measured to be 4.4° for 
a detuning of A 2 = 27r • 850 MHz. This translates into an optical depth of d « 168. 

Having defined the optical depth, we can now estimate the number of classical photons that need to be filtered from 
the quantum photon. The field at the detector (see Fig. lb in the article) is described by the operator d in Eq. (8) in 
the article. Assuming a length, t-mt, of the write pulse, we estimate the average number of quantum photons, A(^uant 
at the detector as 


A^quant — ■ tint) — ^ ^ I )e 


16 


(52) 


where we have used that {\9j\^)e is independent of time as shown in Eq. (14) in the article. Note that {\9j\^)e oc 
\g\^ |r2|^ and we estimate the number of classical photons contained in the write pulse as A^cias ~ tintKi/(4 Ig]^) = 
tintKi/(4/3 \g\^), where (3 = \gg\^ / is the ratio between the Clebsch-Gordan coefficients (/x) of the transitions 


characterized by g and g (see Fig. 8b). From Eq. (52), we then get 


Nclas 


N, 


quant 


^ \ 9 rmiP 4 


(53) 


The number of classical photons that needs to be filtered is finally estimated by setting Nquant = 1- Using Eqs. ( |47| )- 
(48), we can express Adas in terms of the optical depth and the finesse of the cell-cavity, defined as = 27r/(rKi), 
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where r is the cavity roundtrip time. Furthermore, we assume that {\0j'f')e 
photons can be estimated as 


87r/3|L2A2 1 

3/3^27(71 + 72) ’ 


\{dj)e'f' such that the number of classical 


(54) 


where we have expanded the expression for \{0j\ 

the Clebsch-Gordon coefficients of the transitions characterized by 51,52 and 5 in Fig. | 8 b 

7.4-10" 


1 ^ in the limit of large detuning. P 2 = is the ratio between 


For the experimental Cs- 
With d = 168 and T = 100 this gives 


cells and a detuning of A = 27r • 898 MHz, we find that fVcias 
Ncias = 4.4 • 10®. Since the quantum and classical field differ both in polarisation and frequency, this level of filtering 
is expected to be easily achieved using using a combination of both polarisation filtering and the filter-cavity. 


IX. READOUT 


In this section, we give the expressions for the second order correction to the readout efficiency and present the 
details of the simulations. For the second order term reread ,2 we find 


2|I?| 


5read.2 = 2«;iReal( / dt I dt' I 

\ Jo Jo Jo 2 \V\ \ 




Vnb^ 
Vv I 




A-C+Vv\ 


2VV J 

|—.A-l-C -f 'J^D\ 
2VV ' I 


+e 


+e 


'2 2 










.A-C-Vv\ 

^ 2VV ) 

-A+c-Vv) 


2y/v 


(55) 


Here we have, once again, neglected the contributions from the fluctuations contained in 5A{t) and Bq since they are 

^( 2 ) 
•‘cell 


suppressed by a factor of at least dT/N compared to the terms above. In deriving Eq. ([^, we have used that 


consists of sums of the form 


I 

N 


N-l N N 


E E E 


(56) 


1=1 i=ii'=i 


where Xj could e.g. denote Bj. For 5 read, 2 , we calculate (d^eii^ceii) average of Eq. (56) is approximatively 

N{dXj{t')SXj{t''))e because {6Xj{t')6Xji{t'')) = 0, if j ^ j', since the motion of different atoms are uncorrelated and 
we have assumed that A — 1 « iV. All correlations appearing in Eq. (55) are thus single atom correlations and the 


index j is kept to indicate this. The correlations contained in 5 read ,2 can be treated analytically, in a similar fashion 
as the correlations in {\dj{t)\^) for the write process, but we have instead simulated the correlations numerically for 
the previously mentioned Cs-cells. 


A. Numerical simulation 

The simulations are performed in the same way as for the write process. An extra difficulty is, however, that we 
consider the coupling between the light fields and the extra levels in ®®®Cs. We assume that the readout process has 
the level structure shown in Fig. 
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FIG. 9. Schematic view of the readout realized in the hyperfine levels of 6^5'i/2 and 6 ^P 3/2 in We imagine the single 

exitation to be stored in state \F = 3,mF = 3) while the macroscopically populated state is \F = 4,ttif = 4). Note that 
the classical drive also couples \F = 4,= 4) to \F' — 4:,mpi = 4) and \F' = 5,mpi = 4), which can pump atoms out of 
\F = 4, rriF = 4). These couplings are however sufficiently suppressed by the large splitting of 2tt ■ 9.2 GHz between the ground 
states. 


The couplings to the extra levels result in extra coupling terms in the expressions for A,Bj 


and Cj, 


which we 


include, but the expression for reread ,2 is still the same as given in Eq. (55). Note, however, that a cavity detuning of 
the quantum field (appearing in the expression for A) on the order of Ki is needed to compensate the phases resulting 
from some of these additional couplings. The starting point of our numerical simulations is therefore Eq. (55), where 
we can change the order of integration and introduce the variables u = t' + 1" and s = t' — t" since the correlations 
only depend on the time difference \t' — Performing the integrals over t and u analytically, allows us to write 


^read,2 — 


(/il('rread,s)(dS,(5;B)e(s) + ^2 (r^ead, s) (<5^, e (s) 


Th3 (Tread) s) (dC, e (s) 4“ ^4 (H'ead, s) (JC, (5C) g (s))ds, 


(57) 


where hi(Ti.eadi s), h 2 (Ti.ead) s), h 3 (Tread) s) and /i 4 (Treadj s) are functions of s and Treadj which are obtained from the 
integration over t and u. We have once again introduced the short notation for the correlations {5Bj{t')5Cj{t"))e = 
{SB, SC)e{s). Note that {SB, SC}e(s) —>■ 0 for s —^ oo similar to the situation in the write process, i.e., the coupling of 
an atom at time t is uncorrelated from its initial coupling if t is large. We can therefore introduce a cutoff Smax in 
the integral in Eq. ( |h^ such that we can evaluate reread ,2 for an arbitrary length of the readout pulse Tread without 
additional numerical difficulty. We then numerically evaluate reread ,2 from Eq. (57) by simulating the decay of the 
correlations similar to the simulation of the write process. From the numerical simulation, we find that the term 
{SCjSCj)e dominates 77read,2- This term describes loss of the excitation due to spontaneous emission to modes not 
confined by the cavity. 


X. ERRORS 


Here we give the detailed expression for the probability (pi) to read out incoherent photons to first order, 
incoherent photons are being readout from assymetric modes described by the operators 


The 
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Si = 


.7=1 


„2ipU-pl/NAj) 


mi ! 


(58) 
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where Z = {1, 2,..., TV — 1}. Note that the symmetric Dicke mode is described by Sq. From the pertubative expansion 
of Oceii we find that the assymetric modes gives a first order contribution of 


j(i) 

^cell 


(t) = ^/€ / dt' 


t g5(^+C)(t-t') 2 


7V-1 N 


2^/V VN 




-2iTr{j-l)l/N j 


_g 2 


1 = 1 3 = 1 

{{A-C) + 2B5Cj{t')) Si 


+ y/V6Bj{t' 

The probability to read out the incoherent photons is then 


Pi = 




dt f dt' f (t'))^aL\(f"))- 

JO JO JO 


(59) 


(60) 


We can numerically evaluate the correlations contained in Eq. (60) in a similar fashion as the correlations in reread ,2 j 
i.e., we simulate the experimental Cs-cells. Thereby we get the results presented in Fig. 3b in the article. 


A. DLCZ error analysis 


In order to fully characterize the performance of a DLCZ repeater based on the room temperature cells, we inves¬ 
tigate how the errors from incoherent photons propagate in the repeater. First we consider the state being produced 
in the entanglement generation step (see above). We define the following parameters: 

• ^write : The write efficiency, which basically is the probability to have a Dicke state in the ensemble conditioned 
on a quantum photon being emitted. If we are not in the right state we are dominated by the intitial state and 
thus with probability 1 — r/^rite) we assume the atomic state to be |00... 0). 

• Pe '■ The excitation probability, which depends on the driving strength. 

• Pd : The dark count probability of a detector. We estimate this as Pd ~ rdarkUnt where r^ark is the dark count 
rate and tint is the length of the driving pulse. 

• rjd : The detection efficiency, which is determined by the total efficiency of the detector, the outcoupling losses 
and the transmission losses from the cavity to the detector. 

To second order in the excitation probability, the state of the two ensembles following a single click in a detector 
at the central station is described by the density matrix 


Pe(l - " -Pd)" + 2pl{l - r;write)"(l " 77d)"-Pd(l " -Pd) 

-1-2(1 -pe)"P’d(l - Pd) + 4pe(l - 77write)(l - Pe)A “ Vd)PdA “ -Pd) 

+4p2(l _ r;write)"(??d - r?d)(l - -Pd)" + 4p2(I - ?7write)"(l " r?d)"-Pd(l - -Pd)] | 00 )( 00 | 

+ PlvLitei‘^Vd - 2?7d)(l - Pd)" + 2plrjl^.^^{l - ?7d)"Pd(l - Pd) |11)(11| 

T Pe^write(l ^write)(2?7d 27yd)(i- Pd) T 2pg(1 ^write)(l Pd) Pd(l Pd) 

T2pe^write(l Pe)(l ^d)Pd(l Pd) “t” ^PePvjriteA ^write)(^d ^d)(i^ Pd) 
+4p2r?write(l - d7write)(l - d/d)"Pd(l - Pd)] (| 01 )( 01 | + | 10 )( 10 |) 

-f 


^plvlriteiVd - vDa - Pd)" + 2plpl,,,,il - r,d)"Pd(l - Pd) (| 20 )( 20 | + | 02 )( 02 |) 
+ 2ped?write(l -Pe)?7d(l “ Pd)" I’®') ('®'I 


1 


(61) 


where we have defined |0) = |00...0), |1) = |Dicke), |2) = ^ |l)*|l)i(0|i(0|j |00 ■ • • 0) and |’®') = ^(|01)-f 

|10)). Af = Psuccess is a normalization constant, which gives the success probability of the operation. Note that we 
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have assumed number resolving detectors. It is seen from Eq. (611 that we can write p. 
co(|01)(01 


success 


= ao|’®')('i'|+&o|00)(00| 


|10)(10|) +do|ll)(ll| + eo(|02)(02| + |20)(20^ 

We now consider the entanglement swapping of two states of the form Psuccess using the setup shown in Fig. In 
the swap, an ensemble from each entangled pair is read out and the corresponding photons are combined on a balanced 
beam splitter. With probability Cg, we are swapping two states of the form |’®') and the state after a successful swap 
is 


Pal 


1 r 
4 i 

1 


2r?dr(l - r?dr)(l - Pdrf + 2(1 - VdrfPdril - Pdr) |00)(00 


2Pd,(l-Pd,)J|ll)(ll| 

2(1 - 77dr)^dr(l - Pdr)] (|01)(01| + |10)(10|) 


+ 2^dr(l-Pdrri'P)('&| 


1 

77^’ 


(62) 


where we have introduced the total readout detection efficiency pdr, which is determined by the readout efficiency, 
the outcoupling losses and the efficiency of the detectors. In contrast to the detection efficiency in the entanglement 
generation, the readout detection efficiency does not include fiber losses since the entanglement swap is a local 
process. Af' is a normalization constant. Furthermore, we have defined the readout dark count rate Pdr, which 
contains the probability of reading out incoherent photons from the ensembles and the detector dark counts, i.e. 
Pdr = J'dark'Tread + VdrPi- Here TdarkTread Is the dark count rate of the detectors and pi is the probability to emit 
incoherent photons. As previously described pi is mainly determined by the inefficiency of the optical pumping in the 
initialization of the ensembles and the memory time of the ensemble. We will neglect the limited memory time and 
simply assume a fixed value of pi from the inefficiency of the optical pumping. From Eq. (62), we can express Pg 2 as 


1 

Po'l ~ Jp 


ia|00)(00| + i(/3 - pdr{l - Pdr)")(|01)(01| + |10)(10| 


+ iA|ll)(ll| + ir?dr(l-Pdr)^|^)(’®'| 


(63) 


we find that the output state can be written as 


with a, ,5, and A given by Eq. (62). Considering all the combinations from swapping two states of the form Psuccess, 
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where we have defined 
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FIG. 10. Sketch of the postselection procedure after a single step of entanglement swapping (middle stations). The entanglement 
swaps creates entanglement between ensembles 1,2 and 3,4, respectively. Ensembles 1,2,3 and 4 are then read out and we 
condition on a single click at both stations. The phases </> and 9 are ideally equal. 


and Af" = Jljwap,i is a normalization constant, which gives the success probability of the swap operation. It is seen 
that we can write 


Pswap,i = + 6i|00)(00| +ci(|01)(01| + |10) (10|) + dp |11) (11| + ei(|02) (02| + |20)(20|) 

+/i(|12)(12| + |21)(21|) +5 i|22)(22| (68) 


with constants oi, 6i, ci, di, ei, /i and gi determined by Eq. (64). 


The vacuum part of the swapped state grows exponentially with the number of swaps in the DLCZ protoc ol [3 and 
it is therefore necessary to perform a final postselection where two entangled states are combined (see Fig. 10). We 


assume two parties named Alice and Bob who share two entangled pairs such that Alice has ensemble 1 and 3 and 
Bob has ensemble 2 and 4. Ensemble 1 and 2 are entangled and so is ensemble 3 and 4. The intial state is then 


dintiai = + 6|00)(00| +c(|01)(01| + |10) (10|) + d|ll) (11| 

+e(|02)(02| + |20)(20|) + /(|12)(12| + |21)(21|) + ff|22)(22|) 

/ 1,2 

(g)(o|’i')(4'| +6|00)(00| +c(|01)(01| + |10)(10|) + d|ll)(ll| 

+e(|02)(02| + |20)(20|) + /(|12)(12| + |21)(21|) + ff|22)(22|) (69) 

/ 3,4 

All ensembles are now readout in, e.g., a cryptography scheme and a success is conditioned on both Alice and Bob 
recording a single click. The total success probability is 

^2 

Pps = "^(<aps/3ps + 7ps) 4“ ^ '^ps 4“ 2c Ops/Sps + d /3pg + 2c 7pg + 2abaps^ps 

4“2(zc(ops/4ps T 7ps) T 2nd/4ps7ps 4“ 46co;psyps -t- 4cd/3ps7ps 4“ 26dypg 

4“2tte(ops/4ps 7ps/4ps) 4“ 2cif Pps) 4“ 2ct^/4ps/4ps -)- AbccXpsPps 

4“46/yps/?ps + 26g/4pg + 4ce(yps^ps 4- cxps^ps) 4- 4c/(/3pg + 7ps,dps) 4“ 4c5/3ps/3ps + ddeyps^ps 

4“4d/,dps,0ps + 2dgPpg + 2e (cipsTps 4- /4ps) + 4eg/3ps7ps 4- 4:f gPps^ps 4- g 7ps 

4-4/e(7ps7ps + Ppsldps) + 2/^(7ps^ps 4- (70) 


Ops = 2Pdr(l — d^dr) (71) 

/3pe = 277dr(l - r7dr)(l " ^dr)" + 2(1 - 77dr)"Pdr(l " ^dr) (72) 

7ps = ddr(l - d^dr)^ 4“ 2(1 - ?7dr)d^dr(l “ ^dr) (73) 

Pps = 3?7dr(l - r?dr)Xl - Pdr)^ + 2(1 ” ddr)^Pdr(l " Pdr) (74) 

= 477dr(l - r?dr)Xl - Pdrf + 2(1 - ?7dr)^Pdr(l " Pdr) ■ (75) 


where we have defined 










22 


The postselected fidelity of the state is 


p — 
ps — 


Thread “ ^’dr)'‘(l + COs((/) - 9)) + ^(7^, - - PdvY) + + aCJ^ 




ps 


(76) 


where the phases (f>, 9 can be different due to variations in the path lengths of the photons being read out. In the 
ideal case, we have <j) = 9 for which the fidelity is maximal. 

The rate of a DLCZ repeater based on the room temperature cells can be estimated as [7] 


r 



n+1 




0 wap, 1 wap, 2 


P P 

• ■ ^ swapjTi.* ps 


2 "c 

~17 


(77) 


where n is the number of swap levels in the repeater and L is the distance. Note, however, that the scalability of the 
room temperature cells enables spatially multiplexing, which both increases the rate of the repeater and decreases 
the necessary memory time of the ensembles [7] . Assuming that 2M ensembles are used at each repeater station, the 
rate will increase by a factor of M. We have considered a basic repeater segment consisting of only a single swap 
without multiplexing in order to estimate the rate and fidelity of a distributed pair. We assume that the distance to 
distribute entanglement over is L = 80 km and that the losses in the fibers are given by the attenuation length at 
telecom wavelengths, which is ~20 km. Furthermore, we assume SPD efficiencies of 95% and dark count rates of 1 
Hz. This reflects what is possible with current superconducting detectors mi- We then perform an optimization in 
all the parameters characterizing the cells, e.g. the excitation probability, the write time and the readout time. We 
include experimental imperfections such as outcoupling losses of around 10% and intracavity losses of 2%. As a result 
we find that a pair with fidelity ^ 80% can be distributed with a rate of ~ 0.2 Hz. We have assumed that e « 0.5% 
and have neglected effects from finite memory time of the atoms. Furthermore, we have assumed that (j) = 9 (see 

Fig.ra. 
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